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ION RADIAL DIFFUSION IN AN ELECTROSTATIC IMPULSE MODEL 
FOR STORMTIME RING CURRENT FORMATION 


Margaret W. Chen, Michael Schulz, Lany R. Lyons, and David J. Gomey 
Space and Environment Technology Center, The Aerospace Corporation 


Abstract. Guiding-center simulations of stormtime 
transport of ring-current and radiation-belt ions having 
first adiabatic invariants p £ 15 MeV/G (E ~ 165 keV 
at L ~ 3) are surprisingly well described (typically within 
a factor of i. 4) by the quasilinear theory of radial dif- 
fusion. This holds even for the case of an individual 
model storm characterized by substorm-associated im- 
pulses in the convection electric field, provided that the 
actual spectrum of the electric field is incorporated in the 
quasilinear theory. Correction of the quasilinear diffusion 
coefficient D\ L for drift -resonance broadening (so as to 
define PrV) reduced the typical discrepancy with the dif- 
fusion coefficients deduced from guiding-center sim- 
ulations of representative-particle trajectories to a factor 
~ 3 The typical discrepancy was reduced to a factor 
- 1.4 by averaging Dff , D\ l , and D r A over an en- 
semble of model storms characterized by different (but 
statistically equivalent) sets of substorm-onset times. 

Introduction 

In a guiding-center simulation designed to investigate 
the access of particles to form the stormtime ring current 
[Chen et al ., 1992], hereafter referred to as paper 1, we 
have found that particles having energies E < 110 keV 
could have gained access to L *** 3 essentially by direct 
convection [Lyons and Williams, 1980] from the plasma 
sheet, whereas the inward stormtime transport of parti- 
cles that attain energies E £ 165 keV at L ~ 3 resembles 
radial diff usion [cf. Cornwall , 1968; Lyons and Schulz , 
1989]. We made quantitative comparisons of diffusion co- 
efficients obtained from the simulation with those derived 
from quasilinear diffusion theory [Fdlthammar, 1965] in 
an effort to show that the transport is indeed appropri- 
ately described as radial diffusion. The diffusion coef- 
ficients constructed from simulated particle trajectories 
agreed surprisingly well (typically within a factor of ~ 4) 
with the predictions of quasilinear radial-diffusion the- 
oiy ? even over a model storm whose main-phase duration 
(taken as 3 hr) amounted to only about 3-30 drift periods 
for representative particles of interest. 

In the present work we explore two refinements which 
are expected to produce theoretical diffusion coefficients 
that agree even better with those obtained from sim- 
ulated particle trajectories. The first is a resonance- 
broadening correction [cf. Dungey, 1965] to the quasilin- 
ear diffusion theory of Fdlthammar [1965] for an individ- 
ual model storm. The second is an application of quasi- 
linear theory to the mean diffusion coefficients obtained 
by simulating particle trajectories over an ensemble of 20 
statistically similar model storms. 

Field Model 

We use the magnetic field model that is obtained by 
adding a unif orm southward field AB to the geomagnetic 
dipole field. This superposition leads to the appearance 


of a quasi-magnetopause at the bou ndar y between closed 
and open field lines. The equation of a field line in this 
model is 


[1 + 0.5(r/6) 3 ] _1 (r/f2f:)csc 2 B = constant = L (1) 

where ris the geocentric distance, 0 is the magnetic colat- 
itude, Re is the radius of the Earth, and b = 1.5L Rs = 
12.82 Re is the radius of the equatorial neutral line. This 
value of b which is obtained by mapping the last closed 
field line (denoted L *) to a colatitude of 20° on the Earth, 
corresponds to |AB| = 14.474 nT and L = 8.547. The 
limit b oo (I* -* oo) would correspond to a purely 
dipolar B field. In this study, we consider only equa- 
torially mirroring particles. The equatorial field inten- 
sity Bo is given by Bo = (mf/ 7 * 3 ) — 14.474 nT, where 
fi E = 3.05 x 10 4 nT -R? E is the geomagnetic dipole mo- 
ment. Further details of this field model are given by 
Schulz [1991, pp. 98-110]. 

We assume that the total electric field E = -V$£ is 
from the scalar potential 




( 2 ) 


in which the three separate terms correspond to corota- 
tion (Vn = 90 kV), the Volland-Stem model of quiescent 
convection (Vo = 50 kV), and the time-dependent en- 
hancement AV(f) associated with the storaitime convec- 
tion, respectively. The time- varying term in the potential 
is assumed to vary as L rather than as L“ because electric 
disturbances are expected to be less well shielded than 
steady-state convection by the inner magnetosphere. 

We model the storm- associated enhancement A V(t) m 
the cross-tail potential drop, 


N 


A V(t) = Y, AVi exp K f * “ ~ *«) 


(3) 


*= l 


where 0(f) is the unit step function (= 1 for t > 0; =0 
For t < 0), as a superposition of almost randomly occur- 
ring impulses that rise sharply and decay exponentially 
with a “lifetime” r = 20 min [cf. Cornwall, 1968]. The 
impulses represent the constituent substorms of a storm. 
The potential drop A Vi associated with any individual 
impulse is chosen randomly from a Gaussian distribution 
with a 200-kV mean and a 50-kV standard deviation. We 
assume that the N start times t, in (3) are randomly dis- 
tributed within a 3-hr time interval corresponding to the 
main phase of a storm, except that we impose a 10-min 
“dead time” (after each impulse onset) during which no 
subsequent impulse can start. This constraint imposes 
a realistic delay between the occurrences of consecutive 
impulse onsets. Without such a dead-time it would be 
possible for the next impulse to start immediately after 
the previous one, and this could lead to the build-up of 
unrealistically large cross-tail potentials. Further details 
of this model storm are given in paper 1. 

We have constructed 100 such random storms so that 
on average there are 9 impulses per storm or 3 sub- 
storms /hr. We have randomly chosen one model storm 
for a detailed case study. Figure 1 shows the variation in 
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cross-tail potential for this prototypical storm. The mean 
enhancement in the cross- tail potential drop for this par- 
ticular storm over the time interval t\ < t < t\ -f 3 hr is 
(AVXO) = ISO kV. 

Simulated Trajectories 

In paper 1 we considered the stormtime transport of 
singly charged ions with various y values by numeri- 
cally tracing their guiding-center trajectories (Burlish- 
Stoer extrapolation method). We began by tracing the 
quiescent trajectories obtained by setting AV(t) = 0 in 
(2). The solid curves in Figure 2 illustrate steady-state 
trajectories of equatorially mirroring ions for a particular 
H value (30 MeV/G). This corresponds to an energy of 
335 keV at a geocentric radial distance of r = 3 Re- For 
detailed study we chose the quiescent drift shell that in- 
tersects the dusk meridian at R = t/Re = 3 for each of 
12 y values (15 MeV/G 200 MeV/G). We selected 

this drift shell because it is representative of where par- 
ticles need to be transported in order to form the storm- 
time ring current. We computed the quiescent ionic drift 
period T3, which is 0.73 hr for y = 30 MeV/G. We spaced 
12 representative ions equally in time on the drift path 



Fig. 2. Solid curves: quiet-time equatorial trajectories of 
singly charged ions having y = 30 MeV/G. Outer dashed 
circle: neutral line at r = b. Small “filled” circles: the 12 
representative ions 1 “final” positions on the steady- state 
drift path of interest. Stormtime trajectories computed 
in our time-reversed simulation span the stippled area. 


of interest before starting the simulation. To investigate 
particle access to the drift shell of interest, we ran the 
simulation backward in time. Thus, the time- reversed 
trajectories (which span the stippled area in Figure 2) 
indicate where the particles must have been prior to the 
storm in order to have reached the final drift shell (L ~ 3) 
of interest. Indeed, we found a spread among the initial L 
values of the representative ions. Such particle transport 
among closed drift paths resembles radial diffusion and 
was typical for ions with y ~ 15 MeV /G at L ** 3. We 
also ran simulations forward in time, so as to follow the 
dispersal of ions from a common initial drift path. Not 
surprisingly, we found a qualitatively similar transport, 
although the particles followed in the time-forward sim- 
ulations were implicitly different from those followed in 
the time-reversed simulations. 

Diffusion and Quasi- Diffusion Coefficients 

Since the quiescent drift shells are not circular (cf. Fig- 
ure 2), we found it convenient to label drift shells (and 
thus to express the diffusion coefficients) in terms of the 
dimensionless third adiabatic invariant defined by Rotd- 
trer [1970, p. 107] as 

L = \2* l i E \- \tef w>)\' K) 

where $5 is the magnetic flux enclosed by that drift shell 
and L(4>) denotes the field-line label at longitude <t> on the 
drift shell. We computed diffusion coefficients for each 
y from the distribution of initial L values of the time- 
reversed trajectories (see Figure 2) by constructing the 
quantity 



where L% and I>f denote the drift-shell labels of the initial 
and fin*! trajectories, respectively, of the 12 representa- 
tive ions and T (= 3 hr) denotes the duration of the main 
phase of the model storm. The quantity ™ is thus a 
measure of the variance among the initial third adiabatic 
invariants of particles situated on the final drift shell of 
interest. The result for each y is plotted as an open cir- 
cle in Figure 3. We also computed diffusion coefficients 
DVT from the time-forward simulations by interchanging 
the indices i and / in (5), and these results are shown as 
filled circles in Figure 3. The diffusion coefficients DVf 
obtained from simulations run forward and backward in 
time are not very different, 80% of them being within 
24% of each other, although they pertain implicitly to 
different sets of particles. 

We compare the diffusion coefficients obtained from 
the simulated trajectories with the resonant -particle for- 
mulation [F&lth&tniTuiT, 1965] of radial-diffusion theory in 
which the diffusion coefficient is of the form 

where E(u/2ir) is the spectral-density of the (quasi-uni- 
form) equatorial electric field in the inner magnetosphere 
and ft 3 /2x is the particles’ quiescent drift frequency. 
When we substitute the spectral-density function for our 
model storm (see paper 1 for derivation) into (6), we ob- 
tain the quasilinear diffusion coefficient 

ql t 2 L 6 R 2 e A ^ AVj&Vj cos Wj ~ t.)] 

“ ~ l6Tp 2 E (L*) 2 1 + air 2 (7) 
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p y MeV/G 

Fig. 3. Diffusion coefficients Df™ obtained via (5) from 
time-reversed (open circles) and time-forward (filled cir- 
cles) simulations, for comparison with the quasilinear dif- 
fusion coefficient D\ L (dashed curve) given by (7) as an 
implicit function of p and with r (solid curve), which 
represents a correction of the quasilinear diffusion coeffi- 
cient for resonance-broadening effects. 


in which correlations between the impulses lead to cross 
terms (J ^ i). If we neglect the cross terms, we recover 
essentially the standard diffusion coefficient of Cornwall 
[1968]. 

The quasilinear diffusion coefficient at £ ft* 3, repre- 
sented by the dashed curve in Figure 3, is not a veiy 
smooth function of p. This is because the impulse onsets 
(t) associated with any individual storm modeled by (3) 
occur at specific (although randomly determined) times 
U so that the corresponding spectral density E{oj/2r:) is 
not a very smooth function of frequency. The diffusion 
coefficients obtained from the simulations can be seen to 
agree surprisingly well with the quasilinear results despite 
the strong variability of D q LL with p. Typically (i.e., in 
> 80% of the available comparisons) we find agreement 
within a factor of 4. Quasilinear theory even accounts 
for the p values (e.g., p = 75 and SO MeV/G) for which 
the diffusion coefficients computed from the simulation 
are especially small. 

However, the diffusion coefficients obtained from the 
simulated trajectories do not show quite as much vari- 
ability with p as quasilinear theory predicts. This is 
not surprising since quasilinear theory postulates a per- 
fectly sharp resonance at the quiescent drift frequency, 
whereas the simulated transport leads to an eventual 
spread among the quiescent drift frequencies of the rep- 
resentative ions for each p. A rough estimate for the 
anticipated spread in ^lz/2^ is 

^l2**{D LL TI2* 2 ) ll2 \{dSlzldL)r\ ( 8 ) 


Since the frequency bandwidth given by (8) depends on 
Dll, we have iterated between (9) and (8) until satisfac- 
tory convergence to the desired solution (called D^r) is 
achieved. The result is plotted as a solid curve in Fig- 
ure 3. We find that the inclusion of nonlinear resonance- 
broadening effects reduces the discrepancy between quasi 
linear theory with D*™ through a smoothing of the sharp 
relative maxima and minima (with respect to p). Correc- 
tions were typically — 10 - 30% at the relative maxima 
but were as much as 60% at the relative minima (e.g., 
near p — 184 MeV /G). It is not clear whether the remain- 
ing, discrepancy (typically a factor £ 3) between 
Dll ’ * s intrinsic to this kind of study or attributable to 
the neglect of other nonlinear and/or quasilinear effects. 
A preliminary test suggests that the replacement of L 
in (9) by its transport- averaged value (L 6 ) would also 
be a relatively unimportant nonlinear correction. How- 
ever, a quasi-linear correction for the variations of L{4>) 
and <)> along a quiescent drift trajectory would introduce 
harmonic resonances (a? = mf^) in (9) and dim in ish the 
weight of the fundamental drift resonance (w = &). We 
hope to investigate this correction in a future study. 

Diffusion Averaged over an Ensemble of Storms 

It could be argued that quasilinear theory is more ap- 
propriately applied to an ensemble of model storms than 
to an individual storm. We tested this hypothesis by ran- 
domly choosing 20 different storms from among the 100 
storms that we had constructed. We averaged the diffu- 
sion coefficients obtained from the simulations, standard 
quasilinear theory, and the resonance-broadened quasilin- 
ear theory over the 20 storms. The results are shown in 
Figure 4. The ensemble-averaged quasilinear diffusion co- 
efficient D^l i dashed curve) and its resonance- broadened 
counterpart D r A (solid curve) are considerably smoother 
than Dll and respectively, for an individual storm. 
The ensemble- averaged diffusion coefficients D^™ from 
the time-reversed and time-forward simulation (open cir- 
cles and filled circles, respectively) generally agree much 



p. y MeV/G 


since the mean-square spread in L is 2 Dll'?. An esti- 
mate for the diffusion coefficient (corrected for resonance- 
broadening effects) is thus 



Fig. 4. Ensemble- averaged diffusion coefficients D^T 
obtained from the time-reversed (open circles) and time- 
forward (filled circles), simulations, for comparison with 
quasilinear theory ( Dll ^ function of p , dashed curve) 
and with D T , K (solid curve), which represents ensemble- 
averaged diffusion coefficient corrected for resonance 
broadening effects. 
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better with the theoretical diffusion coefficients in Fig- 
ure 4 than does Df™ with the theoretical diffusion co- 
efficients for pi individual storm in Figure 3. Typically 
&V™ in Figure 4_agree within a factor of only 

1.4. Agreement between Df™ and D r £ L is even better 
(typically within a factor of 1.3). 

Discussion and Summary 

In this work we have considered two refinements of 
qua si li n ear radial diffusion theory for the purpose of in- 
terpreting guiding-center simulation results. We found in 
paper 1 that the quasilinear diffusion theory of Faith am- 
mar [1965] already provides a surprisingly good descrip- 
tion of charged-particle transport at first invariants /i > 
15 MeV/G (E £ 165 keV) around L — 3. This holds even 
for the case of an individual model storm characterized 
by substorm-associated impulses in the convection elec- 
tric field provided that the actual spectrum of the electric 
field is incorporated in the quasilinear theory. 

Here we have found that a self-consistent correction 
of the quasilinear diffusion coefficient Dvr for resonance- 
broadening effects accounts even better to r the diffusion 
coefficients Du obtained from a simulation of particle 
trajectories (reducing the typical discrepancy from a fac- 
tor ~ 4 to a factor ~ 3). Resonance broadening tends 
to reduce the variability of Dll with /x by requiring that 
the electric spectral density E(uj/ 2x) be averaged over 
an appropriate bandwidth Au;/27r surrounding the drift 
frequency fl3/27r. The resonance-broadening correction 
turns out [Schulz, 1975] to be unimportant ( £ 1 %) when 
applied to the standard model [Cornwall, 1968] of ra- 
dii diffusion in which the spectral density is regarded 
as a smoothly varying function of frequency. However, 
the spectral-density for any one of our model storms 
is not smooth at all. Indeed, it typically shows varia- 
tions of at least 2 orders of magnitude over bandwidths 
Au?/2tt ~ 1 hr” 1 , whereas resonance-broadening theory 
typically requires E(w/2n) to be averaged over a band- 
width Ao>*/27t ~ 0.25 hr” 1 . The resonance-broadened 
correction is thus quite significant for our typical model 
storm. 

Even better agreement between D T £ L and DVT is found 
when both are averaged over an ensemble ojT20 model 
storms.. Indeed, the typical discrepancy between D r A 
and (and between M- and Dfi) is only a factor 
~ 1.4, mainly because D q LL and D T ^ are much smoother 
functions of /x for an ensemble of storms than are d ll 
and Djr r for an individual storm. The improved agree- 
ment achieved by ensemble-averaging is consistent with 
the usual interpretation of quasilinear theory (with or 
without the correction of resonance broadening) as a sta- 
tistical description of charged-particle transport. What 
is surprising is that quasilinear theory provides a reason- 
ably good description of the particle transport even for 
an individual storm. 
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After having -made ^he ensemble- average comparison 
between DY™ and D^^, we learned of a recent similar 
study in which P. Riley and R. A. Wolf (preprint, 1991) 
have likewise found better agreement between diffusion 
coefficients obtained from simulated guiding-center tra- 
jectories with quasilinear diffusion theory by averaging 
over an ense m ble of model storms than by making the 
comparison for an individual (in their case, actual) storm. 
However, they restricted their attention to particle ener- 
gies £ 130 keV at L = 3, whereas we have found the 
transport to be quasi-diffusive only for particles that at- 
tain energies > 165 keV at L = 3. 
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